Improved all-order results for the one-loop QED correction to the hyperfine structure in light H-like 

atoms 
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A calculation of the one-loop self-energy and vacuum-polarization corrections to the hyperfine splitting of the 
Is and 2s states in light H-like ions is carried out to all orders in the parameter Za. Using the known values for 
the Za-expansion coefficients, the numerical data obtained are extrapolated from Z — 5 and higher to Z — 0, 
1, and 2, with the resulting accuracy being significantly better than in previous evaluations. Our calculation 
shifts the theoretical value of the normalized difference of the Is and 2s hyperfine- structure intervals in 3 He + 
by 0.056 kHz and improves its accuracy. 
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Introduction 

Hyperfine splitting of the ground state in light H-like sys- 
tems, such as hydrogen, deuterium, tritium, and helium-3 ion, 
has long been known experimentally with extremely high pre- 
cision. The present-day theory of the ground-state hyperfine 
structure (hfs) is still far behind the experiment, due to a rela- 
tively large contribution of the nuclear-structure effects, which 
cannot be accurately calculated at present. One of the possi- 
bilities to overcome this difficulty 1 1| is to study the normal- 
ized difference A21 = 8 V2s — v\ s , where v\ s and v% s are the 
Is and the 2s hfs interval, respectively. A large class of cor- 
rections to vis and v^g (among them, all lowest-order nuclear 
effects) are proportional to the nonrelativistic electron density 
at the position of the nucleus (r — 0) and, therefore, do not 
contribute to the difference A21. Consequently, the theoret- 
ical study of this difference can be performed up to a much 
higher accuracy than that of v\ s and v^ s separately. 

The experimental value of the difference A21 is obtained 
by combining results of two independent measurements of v\ s 
and V2s and is known less precisely than v\ s . The best accu- 
racy is obtained for the helium-3 ion in a combination of two 
relatively old results 001, 



A 2 i( 3 He + ) = 1 189.979 (71) kHz . 



(1) 



Recent progress was achieved in the measurement of v^s in 
hydrogen |4] and deuterium |5], which significantly improved 
the corresponding experimental values for the difference A21. 

Theory of hfs and, specifically, of the difference A 2 i in 
light H-like atoms has recently been examined in detail in 
Ref. 1 6]. It is demonstrated that one of the major uncertain- 
ties in the theoretical prediction of A21 ( 3 He + ) stems from the 
one-loop self-energy correction. The self-energy correction is 
also responsible for a significant part of the theoretical uncer- 
tainty for the ground-state hyperfine splitting in muonium |7|. 

The goal of the present investigation is to improve the nu- 
merical accuracy of the one-loop QED correction for the Is 
and 2s states in light H-like atoms. Our consideration will 
be carried out to all orders in the parameter Za (Z is the 



nuclear charge number and a is the fine-structure constant). 
All-order calculations of the self-energy hfs correction in H- 
like ions have been previously performed by numerous au- 
thors © S IIH E3 EH E3I] . Different evaluations are 
generally in good agreement (except for the first two calcula- 
tions, see a discussion in Ref. 1 14]). For high- and middle-Z 
ions, results for this correction can be presently considered 
as well established at the level of the experimental interest. 
In the low-Z region, however, the experimental accuracy is 
much higher and technical problems encountered in all-order 
calculations are more demanding than for higher- Z ions. It 
would be clearly preferable to perform a direct all-order cal- 
culation for Z = 1 and 2 with an accuracy significantly higher 
than the one obtained from the Za expansion, as it was done 
for the Lamb shift flrjll . However, such a project has not been 
realized yet. Blundell et al. IYM obtained the higher-order (in 
Za) contribution to the Is self-energy correction for Z = 1 
by extrapolating their numerical results for higher Z. Similar 
procedure was employed in the investigation by two of us iTfl 
for the self-energy correction to the Is and 2s hfs intervals for 
Z = 1 and 2. 

The vacuum-polarization hfs correction was evaluated to all 
orders in Za in Refs. HE3II1 (without the magnetic-loop 
Wichmann-Kroll correction) and in Refs. I13l 1 1 9fl (complete 
calculations). However, the above studies were mainly con- 
cerned with high- and middle-Z ions, so that little information 
was provided about the behavior of the Wichmann-Kroll part 
of the vacuum-polarization correction in the low-Z region. 

In the present work, we perform a calculation of the one- 
loop self-energy and vacuum-polarization corrections to the 
Is and 2s hfs intervals in H-like ions. The paper is orga- 
nized as follows. In Sec.|I]we evaluate the self-energy correc- 
tion by employing the additional-subtraction scheme lEoll that 
improves the convergence properties of the resulting partial- 
wave expansion. In this way, we significantly increase the 
accuracy of the numerical results for Z > 5 as compared to 
the previous evaluations. The vacuum-polarization correction 
is evaluated in Sec. [n] The higher-order one-loop QED con- 
tribution is inferred from our all-order results in Sec. IHIIbv 
subtracting the known terms of the Za expansion. Finally, we 
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discuss the experimental consequences of our calculation. 

Relativistic units (h = c = m = 1) are used throughout the 
paper. 



L SELF-ENERGY CORRECTION 

In this section we describe the evaluation of the self-energy 
hfs correction without any expansion in Za. We start with 
general formulas for the self-energy correction in the presence 
of an additional perturbing potential 8V . To the first order 
in SV, the self-energy correction is given by the sum of the 
irreducible, the reducible, and the vertex correction l2lll . 



AE SE = AAr + AE Icd + AE V 



(2) 



The irreducible part arises through a perturbation of the wave 
function, 



AE; r 



<<5a|7 £(£a)|a) + <fl|7°£(e a )|5a}, 



(3) 



where £ — £ — 8m, 8m is the one-loop mass counterterm, £ 
is the one-loop self-energy function, 

/oo 
-oo 

xG(s-LJ,x 1 ,x 2 )a u D^(w,x- L2 ) ,(4) 

G is the Dirac Coulomb Green function G(e) = [e — TL{1 — 
iO)] -1 , TL is the Dirac Coulomb Hamiltonian, D^ v is the pho- 
ton propagator, a M = (1, a), and X12 = X\ — x 2 . The per- 
turbed wave function is given by 



\8a) 



\n){n\SV\a) 



(5) 



The reducible part can be considered as a correction due to the 
first-order perturbation of the binding energy, 



A£ red = 8e a (o|t° |-£>) 

OS 



where 8e a = {a\8V\a). The vertex part is given by 



(6) 



AE V 



2n 



E 



du> 



{ni\8V\ri2) {an 2 \I{uj)\n\a) 



e ni (l - iO)][e a - uj - e„ 2 (l - iO)] 



(7) 



where 7(w) = e 2 a^a v D^(u). 

The self-energy correction to hfs is given by the above for- 
mulas, in which we should assume the perturbing potential to 
have the form of the Fermi-Breit interaction (the nuclear mag- 
netic moment is denoted by fi), 



SV^V hfs (r) 



\e\ a ■ [fj, x r] 



4"7T 



(8) 



and the initial-state wave function \a) to be the wave function 
of the coupled system (electron+nucleus), 



|o> -» \FM F Ij) = 



Mjm, 



c 



FM F 



\IM!)\j a m a ), (9) 



where \IMj) denotes the nuclear wave function, \ j a m a ) is the 
electron wave function, F is the total momentum of the atom, 
and Mp is its projection. Radial integrations over the nuclear 
coordinates can easily be performed already in the general ex- 
pressions. One can show that formulas 0-0 yield correc- 
tions to hfs if we employ the perturbing interaction in the form 



SV(r) 



E F [r x cx] z 



A/i{Zaf 



(10) 



where Ep is the nonrelativistic Fermi energy, and consider 
the initial-state wave function to be the electron wave function 
with the moment projection m a = 1/2, 

\a) = \j a l/2). (11) 



A. Irreducible part 

As follows from Eq. 0, evaluation of the irreducible part 
of the self-energy hfs correction implies a calculation of a 
non-diagonal matrix element of the self-energy function and, 
therefore, is very similar to the evaluation of the first-order 
self-energy correction to the Lamb shift. Since our present 
approach to this problem is somewhat different from the stan- 
dard potential-expansion method, we now give a short de- 
scription of the scheme used for the evaluation of a self-energy 
matrix element. 

Ultraviolet divergencies in the self-energy function are 
traditionally isolated by separating the first two terms in the 
expansion of the bound-electron propagator G in terms of the 
binding potential V, 

G(E,xi,x 2 ) = G^(E, Xl ,x 2 )+G^{E, Xl ,x 2 ) 

+G {2+) (E, Xl ,x 2 ), (12) 

where = [w - Ho(l - iO)]" 1 is the free Dirac Green 
function, G^ is the first-order expansion term 

G {1) (E, Xl ,x 2 ) = J dz G (Q) (E,x 1 ,z)V(z)G (0) (E,z,x 2 ). 

(13) 

and G^ 2+ - ) is the remainder. Representing G in the form dl 2l i 
leads to splitting the matrix element of the self-energy func- 
tion into the zero-potential, one-potential, and many-potential 
parts (see Refs. 1221123112411 for details), 

(o|7 E(e o )|o> = AE zcm + AE onc + A£ many , (14) 

with the mass-counterterm part naturally ascribed to the zero- 
potential term. 
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AE 



many 



Modifications of the standard potential-expansion approach 
introduced in our previous investigation 12011 concern the 
many-potential term, which is given by 

2ia / dw j dxidx 2 D^(lu,Xi 2 ) 

J Clh J 

x ^(xi)a M G^ 2+ )(e - u, x 1 , x 2 ) a v ip a (x 2 ) ■ 

(15) 

The integration contour Clh consists of two parts, the low- 
energy (C/,) and the high-energy (G#) one. The low-energy 
part extends from eo — iO to — iO on the lower bank of the 
branch cut of the photon propagator and from iO to eo + iO on 
the upper bank of the cut. In order to avoid the appearance of 
poles of the electron propagator near the integration contour, 
each part of Cl is bent into the complex plane if the calcula- 
tion is performed for an excited state. The high-energy part of 
the contour is Ch — (eo — ioo, £o — iO] + [eo + iO, eo + ioo). 
The parameter e of the contour is chosen arbitrary from the 
interval Eq £ (e a — e% s , 1 + e a ), where e\ s is the ground-state 
energy. 

The function G< 2+ ) =G- G^> - that enters Eq. (Q3) 
is not known in its closed form at present and, consequently, 
an evaluation of the many-potential term has to be performed 
by expanding G (and, therefore, G^ 2+ ^) into eigenfunctions 
of the Dirac angular momentum with the eigenvalue n. This 
expansion will be referred to as the partial-wave expansion 
in the following. The convergence properties of this expan- 
sion are of crucial importance for the numerical evaluation of 
the self-energy correction. In our previous investigation |20], 
it was demonstrated that the convergence rate of the partial- 
wave expansion could be significantly enhanced by separating 
from AE'many a part that is calculated in a closed form. 

It was shown that in the region x\ ~ x 2 the function G^ 2+ ' 

(2+) 

could be approximated by a simpler function G a , 

G^(E,x 1 ,x 2 )=G^(E + il,x 1 ,x 2 ) 

d 



G^(E, 



Xl,X 2 



where 



2Za 



n-^G^(E, Xl ,x 2 ), 
(16) 



(17) 



x\ + x 2 

The above approximation can be obtained from the exact ex- 
pression for G^ 2+ ' by neglecting the commutators [V, G^] 
to all orders. The function G a is expressed in terms of the 
free Green function and can be easily evaluated in a closed 
form. We thus write AE many as a sum of two terms, 



A Tp A Tp sub 

""many — ""many 



■AE 



rcmd 
many ' 



(18) 



The subtraction term AE^^y is obtained from the high- 



energy part of Eq. (1151 by the substitution G' 2 



(24 



The second term AE™%~ is the remainder. The subtraction 
term is evaluated numerically in its closed form (i.e., without 
any partial-wave expansion), whereas the remainder yields a 
rapidly converging partial-wave expansion; for the details of 
the evaluation see Ref . 12011 . 



B. Reducible part 

The reducible part is defined by Eq. (|6j. Using the defi- 
nition of the self-energy function and employing the contour 
Clh for the integration over u>, we write the expression in the 
form 

AE re d = 2 ia 5e a I duj dx\ dx 2 D^(lu, xi 2 ) 

OL v 1p a (x 2 ) ■ 



d 

x ^l{x 1 )a ll —G(e-u,xi,x 2 ) 

OS 



(19) 



Ultraviolet (UV) and infrared (IR) divergences present in this 
expression can be conveniently isolated by separating the 
Green function G into 3 parts, 

G(E) ^ G^{E) + G (a \E)+[G(E)-G^(E)~G^(E)] , 

(20) 

where G^ incorporates the part of the spectral decomposi- 
tion of the bound-electron propagator with e n = e„, 

&>(E,x 1 ,x 2 )-^ E _ £a(i _. oy (21) 

and /i a denotes the momentum projection of the states ip a in 
this expression. The terms on the right-hand-side of Eq. J20l > 
substituted in Eq. dl9l give rise to the splitting of AE^, re- 
spectively, into 3 parts: 



AE, 



red 



In this sum, the term AE^l contains all UV divergences. It is 
calculated in momentum space in a way similar to that for 
the zero-potential part of the first-order self-energy correc- 
tion. UV-divergent terms are covariantly isolated; they dis- 
appear when combined with the free-propagator contribution 
in the vertex part. The term AE^° l contains all IR divergences 
present in the reducible part. (IR divergences of this type are 
sometimes also termed as the reference-state singularities). 
These divergences are regularized by employing the photon 
propagator with a finite photon mass fi. The limit /j, — > 
can be taken when AE^ is combined with the correspond- 
ing contribution from the vertex part. 

The term AE^ ny is finite. Contributions of such type are 
frequently encountered in all-order QED calculations. Usu- 
ally, they are evaluated in coordinate space after expanding 
into an infinite partial-wave series. In the present investiga- 
tion, we modify the standard scheme in order to achieve a bet- 
ter convergence of the partial-wave expansion, analogously to 
that for the irreducible part. We thus separate AE^ 
two parts, 



- AE^ d + AE™* ■ 



(22) 



many 



into 



AE many = ^ 



A^c°d md 



(23) 



The remainder term A£ , I r c ™ ld is obtained from AE' r ™ c j my 
by replacing the standard subtraction [G(E) — G^(E) — 
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(E)] by [G(E) - G^ (E + Cl)- G<» (E)] in the high- 
energy part of the expression. The remaining difference 
[G (0) (£ + CI) - G^(E)] gives rise to the subtraction term 
AE^. More explicitly, the subtraction term is written as 

A£ r s e u d b = 2ia5e a [ dw [ dx l dx 2 D tlu {L0 1 x 12 ) 



-G (0 \e 



d_ r 

de 



■U),Xi,X 2 j 



LJ + Cl, X\,X 2 ) 

a v ip a (x 2 ), (24) 



where CI is given by Eq. < !17i . This expression is calculated in 
its closed form in coordinate space. The calculational formu- 
las are immediately obtained from the corresponding expres- 
sions for the subtraction term for the first-order self-energy 
correction 1 20] . The remainder term is calculated by a partial- 
wave expansion. Due to the additional subtraction in the high- 
energy part, the convergence properties of this partial-wave 
expansion are much better than in the standard approach. 



C. Vertex part 

Rewriting expression (7) for the vertex part of the self- 
energy hfs correction in terms of the bound-electron propa- 
gators, we obtain 

AE VCT = 2ia dio I dx\ dx 2 dx 3 Vj(xi) 

xG(e a - uj,x ll x 2 )SV(x 2 )G{e a -u,x 2 ,x 3 ) 
xa 1/ i) a (x 3 )D^(uj, Xl3 ). (25) 

UV and IR divergences present in this expression can be con- 
veniently isolated by the following separation 

GSVG = G (0) SV G (0) + G (a) SV G (a) 

+ \gSVG— G (0) SV G (0) - G (a) SV G (a) . 

(26) 

This separation, being substituted into Eq. ( I25> . gives rise to 
the following three parts of AE VCI , respectively, 

A£ vcr = AEi°2 + AEg> + AE™r y ■ (27) 

Only the first term in this sum is UV divergent. UV di- 
vergences in AEiel are covariantly isolated by employing a 
momentum-space representation; they disappear when com- 
bined with the corresponding contribution from the reducible 
part, see, e.g., Ref. 0. The second term A£ver is IR diver- 
gent. In order to retain its finite part, we consider it together 
with the corresponding contribution from the reducible part, 

j_ r . i 

2tt 



AEtl + AE, 



(«) 



du> 



E 

Ma'Ma" 



(w - iO) 2 
a'|<5V^|a") (aa \I(ui)\a a) 



-(a\5V\a) {aa'\I(u))\a'a) 



where a' and a" denote the intermediate states with e n = e a 
and with the momentum projection /x / and /Lt a «, respectively. 
The integration over to can be carried out analytically, which 
leads to an explicitly finite result. Sometimes it is more con- 
venient to calculate this contribution directly according to 
Eq. ( 12 8 1 (as long as the contour Clh is employed for the in- 
tegration over oj, this expression is suitable for the numerical 
evaluation). 

The third term AE"^ ny does not contain any divergences 
and is calculated in coordinate space after expanding into a 
partial- wave series. We note that the additional subtraction, 
similar to the one introduced for the reducible part, does not 
improve the convergence properties of the partial-wave expan- 
sion in this case. This is due to the fact that a significant contri- 
bution to the partial-wave expansion terms originates from the 
first-order commutator [V, G^]. In order to achieve a signif- 
icant improvement, one needs to separate the complete con- 
tribution of the vertex with one magnetic and one Coulomb 
interaction. Such contribution was evaluated in a closed form 
for the self-energy correction to the g-factor l25l l26ll using 
the explicit form of the interaction with the constant mag- 
netic field. In the case of the self-energy hfs correction, we 
are presently unable to obtain a closed representation for this 
term. Nevertheless, the partial-wave expansion for AE^ ny 
is converging significantly faster than that for Ai?j" lany and 
AE™^ ny , so that the enhanced convergence achieved for the 
irreducible and reducible parts results finally in a significant 
improvement of the total accuracy of the calculation. 



D. Numerical results 

The self-energy hfs correction is conveniently represented 
in terms of the dimensionless function £>? E defined as 



Ep a 



SE 



D^(Za), 



(29) 



(28) 



where n is the principal quantum number. The results of our 
numerical calculation for the individual contributions of this 
correction for the Is and 2s states and Z = 10 are presented in 
Table H] The calculation was performed for the point nuclear 
model and in the Feynman gauge. In Table [H] we list the final 
results for the self-energy hfs correction for H-like ions with Z 
varying from 5 to 30. A comparison of the results of different 
theoretical evaluations for this correction in the low-Z region 
is given in Table IIIII 

As can be seen from Table IIIII the present calculation im- 
proves the numerical accuracy of the self-energy hfs correc- 
tion by about an order of magnitude for the Is state and even 
more for the 2s state, as compared to our previous calcula- 
tion lfl4ll . This progress is due to the additional subtraction 
scheme employed in the present work for the evaluation of 
the irreducible and reducible parts of the correction. In or- 
der to illustrate the improvement in the convergence proper- 
ties of the partial-wave expansion introduced by this scheme, 
in Figs.[2and|2]we plot (in the decimal logarithmic scale) the 
dependence of the absolute value of the individual terms of 
the partial-wave series on the expansion parameter \k\ within 
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TABLE I: Individual contributions to the self-energy hfs correction for Z = 10, in units of the function D„ defined by Eq. 1361 . 







AE ( T °J d + AfiJW 


AE^l + AfiW 


AB ma„y 




Total 


Is 


-0.263902 


1.896440 


-0.002385 


-1.607827 


-0.185180 


-0.162853(1) 


2 s 


-0.223649 


3.433554 


0.176267 


-1.878883 


-1.649848 


-0.142559(3) 



TABLE II: The self-energy hfs correction for the n = 1 and n = 2 
states of light H-like ions. 



z 


D x (Za) 


D 2 {Za) 


5 


0.174 026 (2) 


0.181940 (2 


6 


0.106 815 (2) 


0.117124 (2 


7 


0.039 476(2) 


0.052 265 (2 


8 


-0.027933(1) 


-0.012 626 (2 


9 


-0.095 379(1) 


-0.077 559 (3 


10 


-0.162 853(1) 


-0.142 559 (3 


12 


-0.297905 (1) 


-0.272 913(1 


15 


-0.501056(1) 


-0.470 078(1 


20 


-0.843 572 (1) 


-0.807153(1 


25 


-1.196 242 (2) 


-1.162 717(2 


30 


-1.566 491 (3) 


-1.547 535 (2 



TABLE III: Comparison of the results of different calculations of 
the self-energy hfs correction for the Is and 2s states in light H-like 
ions, in units of D n (Za). 



Z = 5 


Z = 10 


Z = 20 


Ref. 


Is 0.174 026(2) 


-0.162 853(1) 


-0.843 572(1) 




0.174 028 (20) 


-0.162 860 (20) 


-0.843 588(15) M41 


0.174 05(1) 


-0.162 83(1) 


-0.843 56(1) 


[12] 


0.1741(1) 


-0.162 8(1) 




U3J 


2s 0.181940(2) 


-0.142 559(3) 


-0.807153(1) 




0.18196(10) 


-0.142 51(10) 


-0.807 16 (6) 





the standard potential-expansion approach and within the new 
subtraction scheme. (The parameter k is the Dirac angular- 
momentum eigenvalue of one of the electron propagators in 
the vertex function.) Figs. [^and|2]represent this comparison 
for the irreducible and the reducible part, respectively. 

As a result of the improvement achieved, we were able to 
eliminate completely the uncertainty arising from termination 
of the partial-wave expansion in the irreducible and reducible 
parts. Still, there remains the partial-wave expansion of the 
many-potential vertex term AE™* ny , which has to be termi- 
nated and properly extrapolated to infinity. (In actual calcula- 
tions, the summation was terminated at \k\ = 40). The error 
due to this extrapolation yields one of the main uncertainties 
of our numerical evaluation (the other source of the uncer- 
tainty is the stability of numerical integrations.) Fortunately, 
the partial-wave expansion of A£ , J" 1 any is monotonic and rel- 
atively well-converging (not worse than 1/|k| 3 for all Z > 5), 
and so the uncalculated tail of the expansion can be estimated 
reasonably well. 
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FIG. 1: The absolute magnitude of the individual terms of the 
partial-wave expansion for the irreducible part of the self-energy 
hfs correction for Z = 5, within the standard potential-expansion 
scheme (filled dots) and with the additional subtraction employed in 
the present work (open dots). Plotted are the contributions to the 
function D„(Za) as a function of the absolute value of the relativis- 
tic angular momentum parameter k. A discontinuity of the curve on 
the upper graph around |«| = 19 is due to the change of the sign of 
the contributions to D\ (Za). 



II. VACUUM-POLARIZATION CORRECTION 

The vacuum-polarization hfs correction can be conve- 
niently split into two parts, the so-called electric- and 
magnetic-loop contributions. The electric-loop part originates 
from the diagrams with the hfs interaction attached to the ex- 
ternal electron line, whereas the magnetic-loop one comes 
from the diagram with the hfs interaction attached to the 
vacuum-polarization loop. These contributions are also tra- 
ditionally separated into the Uehling and Wichmann-Kroll 
(WK) parts. The WK part is suppressed by a factor of (Za) 2 
as compared to the Uehling contribution and can often be re- 
garded as a small correction for low-Z ions. 



6 




' i i i i i | i i i i i | i i i i i | i i i i i | i i i i i | i i i i i | i i i i i | i i i i i | r 

5 10 15 20 25 30 35 40 



FIG. 2: The same as in Fig.Q but for the reducible part of the self- 
energy hfs correction. 



The Za expansion of the one-loop vacuum-polarization hfs 
correction reads 



DX P {Za) = {Za)bio + (Zaf 



L b 2 i + b 2 



+ {Zaf Lb 31 +F^ p (Za) 



(30) 



where L = \n[(Za)~ 2 ], the function D^ p is related to 
the energy shift analogously to Eq. i29\ . and the function 
F^ p incorporates all higher-order terms, F^ p {Za) — &30 + 
Za (...)• The first expansion coefficients up to 631 stem from 
the Uehling part of the vacuum-polarization correction; they 
are given by (see, e.g., J(J) 



can be found, e.g., in Refs. lfl3l QjJ. In the case of the 
point nucleus, this contribution was evaluated also analyti- 
cally jg, l27ll - For completeness, we re-calculate it in the 
present work. The corresponding contributions to the higher- 
order remainder F^f p (Za) for the point nuclear model are 
listed in the first and the second columns of Table UVI for the 
Is and 2s states, respectively. The results presented are in 
agreement with the previous calculations of this correction. 

The WK part of the vacuum-polarization hfs correction 
is more difficult to calculate. Especially, this refers to the 
magnetic-loop WK contribution. As outlined in Ref. fl9ll . 
this correction is divergent in the point-dipole approximation 
for the nuclear magnetization distribution. A finite result for 
this correction is obtained if an extended nuclear magnetiza- 
tion distribution is employed. It should be also taken into ac- 
count that the magnetic-loop WK interaction contributes to 
the measured value of the nuclear magnetic moment l28ll . In 
order to prevent double-counting, the corresponding contribu- 
tion should be subtracted from the magnetic-loop WK part of 
the vacuum-polarization hfs correction. Practical calculations 
1H3HI9I1 show that, after such a subtraction, the magnetic-loop 
WK correction depends weakly on details of the nuclear mag- 
netization distribution and has a finite limit in the point-dipole 
approximation. 

The results of our numerical evaluation of the electric- and 
magnetic-loop WK contributions for the Is and 2s states of 
light H-like ions are listed in Table llVI in terms of the higher- 
order remainder F^ p {Za). The electric-loop WK correction 
was calculated for the point nuclear model by employing the 
analytical-approximation formulas for the WK potential from 
Ref. ll29ll . The relative accuracy of this approximation is con- 
sidered by the authors to be not worse than 10" 4 for all Z 
up to Za = 0.95. As an independent test of the accuracy 
of this approximation in the low-Z region, we checked that 
it reproduces well the first two terms of the Za expansion of 
the WK correction to the Lamb shift. The magnetic-loop WK 
correction was calculated for the point-dipole nuclear mag- 
netization model by using a code developed in our previous 
investigation fl^l . A comparison given in Table [TV] demon- 
strates good agreement of our numerical values with the Is 
results of Ref. 1 13] for Z = 10 and 18. 



&io(ns) 



3tt 



(31) 



HI. HIGHER-ORDER ONE-LOOP QED CORRECTION 



621 (ns) 



620 (2s) 
631 (ns) 
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y ' 225 15 



247 
~450 
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(32) 
(33) 
(34) 
(35) 



Higher-order coefficients starting with 630 arise both from the 
Uehling and WK contributions; only their Uehling part is 
presently known lla. 12711 . 

The Uehling part of the one-loop vacuum-polarization hfs 
correction can easily be calculated numerically; some results 



One of the main goals of our investigation is to improve the 
accuracy of the one-loop QED correction for Z = 1 and 2, 
these being the most interesting cases from the experimental 
point of view. The present approach does not employ the Za 
expansion and, therefore, our numerical results do not suffer 
from omission of the higher-order terms, as is the case with 
the perturbative Za-expansion approach. But on the other 
hand, technical problems do not presently allow us to perform 
a direct numerical evaluation for Z = 1 and 2 with a sufficient 
accuracy. In the present work, we employ an indirect method 
used previously in Refs. 1U2I UM . By subtracting the known 
terms of the Za expansion from the all-order results for higher 
values of Z, we identify the higher-order remainder and then 
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TABLE IV: Individual contributions to F^ F (Za) for the Is and 2s states of light H-like ions. 



z 


Uehling 




Electric-loop WK 


Magnetic 


-loop WK 




Is 


2s 


Is 


2 s 


Is 


2s 


1 


7.231 


9.546 


—0.117 


—0.117 






2 


7.337 


9.651 


—0.120 


—0.119 






5 


7.587 


9.901 


—0.128 


—0.125 






10 


7.947 


10.282 


— U.loo 


A 1 OO 

—0.133 


-0.699(2) 


-0.706(2) 












-0.697 a 




12 


8.092 


10.441 


-0.142 


-0.136 


-0.701(2) 


-0.709(2) 


14 


8.240 


10.609 


-0.145 


-0.138 


-0.703(2) 


-0.714(2) 


16 


8.394 


10.788 


-0.149 


-0.141 


-0.705(2) 


-0.718(2) 


18 


8.556 


10.978 


-0.153 


-0.145 


-0.707(2) 


-0.723(2) 








-0.154 a 




-0.706° 




20 


8.725 


11.182 


-0.156 


-0.148 


-0.712(2) 


-0.732(2) 


22 


8.904 


11.401 


-0.160 


-0.151 


-0.717(2) 


-0.740(2) 


24 


9.093 


11.635 


-0.164 


-0.155 


-0.725(2) 


-0.752(2) 



a Ref. 1 13] . 



extrapolate it to Z = 1 and 2. 

First we summarize the results obtained for the one- 
loop self-energy hfs correction within the perturbative Za- 
expansion approach. The corresponding Za expansion reads 



D^(Zu) = a 00 + (Za) oio + {Zdf 



L a 22 + L a 21 + a 20 



+ {Zaf La 31 +F* E (Za) 



(36) 



where L = ln[(Za)~ 2 } and F,f E is the remainder containing 
all higher-order terms, F 7 f E (Za) = a 3 o + Za (...). The 
results presently available for the expansion coefficients are 
(for the references see, e.g., |(J): 



a 00 (ns) 
a w (ns) 
a 22 (ns) 

02l(ls) 
«21 (2s) 

O20 (Is) 

a 20 (2s) 
a 31 (ns) 



In 2- 
2 

3' 



13 

T 



7 

2 ' 



8, „ 37 

-3 ln2+ 72' 
g 

021 (Is) - 3 ln2 " 

17.122 339... , 

O2o(ls) - 5.221233(3). 

5, „ 191 
- In 2 

2 32 



(37) 
(38) 

(39) 
(40) 

(41) 

(42) 
(43) 

(44) 



For the 030 term, there is a preliminary result |30] for the Is 
state, a3 (ls) = —15.9(1.6), and a partial result 15 l3lll for 
the difference A 2 i, a,3o(2s) — <Z3o(ls) = 7.92. 

The higher-order self-energy remainder F„ can be in- 
ferred from our all-order numerical data. The correspond- 
ing results for the function F^ E (Za) and the difference 
F 2 s f(Za) = F 2 SE (Za) - F 1 SE (Za) are plotted on the up- 
per graphs of Figs. [5] and |4j respectively. We note that both 



the F^ E (Za) and F 2 SE (Za) functions have a rapidly vary- 
ing structure in the low-Z region. In order to demonstrate 
this more clearly, we subtract their "linear" part (obtained by 
a global linear fit), with the corresponding plots presented 
on the middle graphs of Figs. and |4] The behavior ob- 
served apparently indicates that the logarithmic term to order 
a(Za) 4 Ep enters with a large coefficient, which complicates 
extrapolation considerably. 

Now we would like to extrapolate our results for the higher- 
order remainder to the lower values of Z, namely Z = 0, 1, 
and 2. For this purpose we employ a procedure similar to the 
one recently described in detail in Ref. 13211 . The extrapolated 
value of a function at Z — zo is obtained in two steps. First 
we apply an (exact) linear fit to each two consecutive points 
from our data set and store the resulting value at Z = zq as a 
function of the average abscissa of the points involved in the 
fit. Then we perform a global parabolic least-squares fit to the 
set of data obtained and take the fitted value at Z — zq as a 
final result. 

We tested this procedure for variation of the logarithmic 
contribution to the next-to-leading order and found it rather 
stable. However, in order to take into account the presence 
of such contribution explicitly, we modify the procedure de- 
scribed above as following. First, we approximate our numer- 
ical data by a function 



f(Z) = c + (Za) [ci ln(Za) + c 2 + (Za) c 3 ] 



with free coefficients 



(45) 



1 by a least- 
squares fit similar to the one described in Ref. [3.3]. Then, 
we use the values obtained for c\ and c 2 in order to define a 
modified higher-order remainder function as 



F° E (Za) = F^(Za) - {Za) [c x \n{Za) 



C2 



(46) 



The numerical results for this function are plotted on the lower 
graphs of Figs. [3] and [4] The function F r f E is much flatter in 
the low-Z region than F n SE and, therefore, is more suitable 
for the extrapolation. We obtain our final results extrapolating 
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-13 




Nuclear charge number Z Nuclear charge number Z 



FIG. 3: The higher-order part of the Is self-energy hfs correction 
_F\ SE as a function of the nuclear charge number (the upper graph); 
-Fi SE with its linear part (obtained by a global linear fit) subtracted 
(the middle graph); F-^ E with its next-to-leading contribution (ob- 
tained by a least-squares fit) subtracted (the lower graph). The nu- 
merical values of the coefficients ci and ci are: ci = 14.83 and 
c 2 = 2.08. 



the function F n by means of the procedure described above. 
The numerical values of the higher-order self-energy remain- 
der obtained in this way are given in the first line of Table IV! 
In the next 3 lines of the table, we present the results of the 
previous evaluations 112, 14, 301. The numerical values ob- 
tained for the function ■F 2 f E in this work fall slightly outside 
the error bars ascribed to our previous results 1 14], as a conse- 
quence of the logarithmic contribution to the next-to-leading 
order being apparently much larger than it was assumed in our 
former work. Our present values for the function F 1 SE are in 
a marginal agreement with the result by Blundell et al. 11211 
and deviate by 1.5 a from the preliminary result by Nio 1 30 ] . 

The Uehling part of the vacuum-polarization hfs correction 
is given in the next line of Table IV! The numerical values are 



FIG. 4: The same as in Fig.^but for the difference = F 2 SE - 
F-±. The numerical values of the coefficients ci and C2 are: ci = 
28.26 and c 2 = 36.85. 

taken from Table|IV]for Z = 1 and 2 and from Refs. BIH 
for Z = Q. The electric-loop WK correction for Z = 1 and 2 
was calculated directly in Sec. [Ill the corresponding numerical 
value for Z = was obtained by a simple extrapolation. Ex- 
trapolation was also employed in order to obtain the results for 
the magnetic-loop WK part of the vacuum-polarization cor- 
rection presented in the table. The error bars specified are ob- 
tained under the assumption that the logarithmic contribution 
to the next-to-leading order enters with a coefficient of about 
2. 

We now turn to the experimental consequences of our cal- 
culation. As demonstrated in Ref. |6], the higher-order self- 
energy correction is one of the major sources of uncertainty of 
the theoretical prediction for the normalized difference of the 
hfs intervals A21 = 8 v% s — v\ s for the 3 He + ion. Our present 
calculation changes the theoretical value of this correction by 
—0.056 kHz (as compared to our former result 1 14]) and im- 
proves its accuracy by a factor of 2. The resulting value of 
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TABLE V: One-loop higher-order QED correction. Acronyms are: "SE" denotes the self-energy contribution, "Ue" - the 
Uehling part, "WK-EL" - the electric-loop WK part, "WK-ML" - the magnetic-loop WK part. 





Fi(Oa) 


Fi(la) 


Fi(2a) 


F 2 i(0a) 


F 2 i(la) 


F 2 i(2a) 


Reference 


jD 




±o.o\o } 


— 14 1 CX\ 




1 .0(1) 










-12.0(2.0) 










[12J 






-14.3(1.1) 


-14.5(7) 




6.5(8) 


6.3(6) 


[14] 




-15.9(1.6) 












[30] 


Ue 


7.06 


7.23 


7.34 


2.32 


2.32 


2.31 




WK-EL 


-0.11 


-0.12 


-0.12 


0.00 


0.00 


0.00 




WK-ML 


-0.69(15) 


-0.69(12) 


-0.69(7) 


0.00 


0.00 


0.00 




Total 


-6.9(4) 


-7.4(3) 


-7.6(3) 


10.7(5) 


9.9(4) 


9.5(3) 





TABLE VI: Normalized difference of the hfs intervals A21 = 8 V2s — 
v\s for the 3 He + ion, in kHz. 



Higher-order QED correction 


-0.594(19) 




A 21, old theory 


-1 190.068 (64) 


[6] 


A21, new theory 


-1 190.124(55) 




A21, experiment 


-1 189.979(71) 





the one-loop QED contribution that incorporates all orders in 
Za starting with the constant term to order a{Za) z Ep for 
the 3 He + ion is given in the first line of Table fVll In the next 
lines of the table, we give the total theoretical value for the 
difference A2i( 3 He + ) taken from Ref. |6], this value modi- 
fied by the present calculation, and the corresponding experi- 
mental result. As can be seen from the table, our calculation 
increases the deviation of the theoretical prediction from the 
experimental value from 0.9 a to 1.6 a. 

It should be noted that our numerical results for all correc- 
tions at Z = 1, except for the one for the magnetic-loop WK 
contribution, can be directly applied to the hyperfine splitting 
in muonium. Our calculation of the magnetic-loop WK cor- 
rection does not hold for muonium since it involves a reg- 
ularization by an extended magnetization distribution of the 
nucleus and a subtraction of the related contribution to the 
measured value of the nuclear magnetic moment. In the case 
of muonium, the nucleus is substituted by a point-like muon 



and the regularization should be performed by a finite mass of 
the muon rather than by a finite size. 



IV. SUMMARY 

In the present investigation, we carried out an all-order (in 
Za) calculation of the one-loop QED correction to the hyper- 
fine splitting of the Is and 2s states in light H-like ions. This 
calculation significantly improved the accuracy of this correc- 
tion, as compared to the previous evaluations. By subtracting 
the known terms of the Za expansion and extrapolating the 
remainder to lower values of Z, we obtained the results for 
the higher-order remainder for Z = 0, 1, and 2. Our calcu- 
lation shifts the theoretical value of the normalized difference 
A21 of the Is and 2s hfs intervals in 3 He + by 0.056 kHz and 
slightly improves its accuracy. 
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